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^^ , Abstract 

t/3 I We deploy algebraic complexity theoretic techniques to construct symmetric deter- 

, ^, ' minantal representations of formulas and weakly skew circuits. Our representations 

produce matrices of much smaller dimensions than those given in the convex geome- 
^^ ' try literature when applied to polynomials having a concise representation (as a sum 

^ I of monomials, or more generally as an arithmetic formula or a weakly skew circuit). 

21 ' These representations are valid in any field of characteristic different from 2. In char- 

QQ I acteristic 2 we are led to an almost complete solution to a question of Biirgisser on 

fT^ ' the VNP-completeness of the partial permanent. In particular, we show that the par- 

tial permanent cannot be VN P-complete in a finite field of characteristic 2 unless the 
(^ • polynomial hierarchy collapses. 

O' 



1 Introduction 

/S ' 1-1 Motivation 

A linear matrix expression (symmetric linear matrix form, affine symmetric matrix pencil) 
is a symmetric matrix with the entries being linear forms in the variables xi, . . . ,rE„ and 
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real number coefficients: 



A{xi, . . . , Xn) = ^0 + xiAi + • • • + XnAn, Ai Symmetric in M ^ . (1) 



A linear matrix inequality (LMI) restricts to those values ^j G M of the Xi such that 
^(?i) • • • 1 ^n) ^ 0, i.e., is positive semidefinite. The set of all such values defines a spectra- 
hedron. 

A real zero polynomial is a polynomial p with real coefficients such that for every x € M"" 
and every /i G C, p{fJ,x) = implies // G M. The Lax conjecture and generalized Lax 
conjecture seek for representations of real zero polynomials /(xi, ..., Xn) as / = det(74) 
with A as in (1) and ^o ^ 0. This is in fact an equivalent formulation of the original Lax 
conjecture which was stated in terms of hyperbolic polynomials (see [Lewis et al. 2005] for 
this equivalence). Furthermore, the matrices are required to have dimension d where d is 
the degree of the polynomial. For n = 2 such representations always exist while a counting 
argument shows that this is impossible for n > 2 [Helton and Vinnikov 2006] (actually, 
[Lewis et al. 2005] give the first proof of the Lax conjecture in its original form based on 
the results of [Helton and Vinnikov 2006]). Two relaxations have been suggested to evade 
this counting argument: At first it was suggested to remove the dimension constraint and 
seek for bigger matrices, and this was further relaxed by seeking for representations of some 
power of the input polynomial. Counterexamples to both relaxations have recently been 
constructed [Branden 2011]. 

Another relaxation is to drop the condition Aq ^ and represent any / as det(^) 
[Helton et al. 2006; Quarez 2008]. However, the purely algebraic construction of [Quarez 
2008] leads to exponential matrix dimensions t. Here we continue the line of work initi- 
ated in [Helton et al. 2006; Quarez 2008] but we proceed differently by symmetrizing the 
complexity theoretic construction by Valiant [1979]. Our construction yields smaller di- 
mensional matrices not only for polynomials represented as sums of monomials but also 
for polynomials represented by formulas and weakly skew circuits [Malod and Portier 2008; 
Kaltofen and Koiran 2008]. Even though in the most general case the bounds we obtained 
are slightly worse than Quarez's [2008], in a lot of interesting cases such as polynomials 
with a polynomial size formula or weakly-skew circuit, or in the case of the permanent, our 
constructions yield much smaller matrices (see Section 4). 

Our constructions are valid for any field of characteristic different from 2. For fields 
of characteristic 2, it can be shown that some polynomials (such as e.g. the polynomial 
xy + z) cannot be represented as determinants of symmetric matrices [Grenet et al. 2011b]. 
Note as a result that the 2-dimensional permanent xw + yz cannot be "symmetrized" over 
characteristic 2 with any dimension. It would be interesting to exactly characterize which 
polynomials admit such a representation in characteristic 2. For the polynomial x -\- y, we 
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where the first matrix is derived from our construction, but the second is valid over any 
commutative ring. It is easily shown that for every polynomial p, its square p'^ admits a 
symmetric determinantal representation in characteristic 2. This is related to a question 
of Biirgisser [2000]: Is the partial permanent VNP-complete over fields of characteristic 2? 
We give an almost complete negative answer to this question. 

Our results give as a by-product an interesting result which was not known to the 
authors' knowledge: Let A be an (n x n) matrix with indeterminate coefficients (ranging 
over a field of characteristic different from 2), then there exists a symmetric matrix B of 
dimensions 0{n^) which entries are the indeterminates from A and constants from the field 
such that det A = det B. This relies on the existence of a size-O(n^) weakly-skew circuit to 
compute the determinant of an {n x n) matrix [Berkowitz 1984; Malod and Portier 2008], 
and this weakly-skew circuit can be represented by a determinant of a symmetric matrix as 
proved in this paper. The dimensions of B can be reduced to 0{n'^) if we replace the weakly 
skew circuits from [Berkowitz 1984; Malod and Portier 2008] by the skew circuits of size 
0{n^) constructed by Mahajan and Vinay [1997]. These authors construct an arithmetic 
branching program for the determinant with 0{n'^) edges, and the arithmetic branching 
program can be evaluated by a skew circuit of size O(n^). After learning of our result, Meena 
Mahajan and Prajakta Nimbhorkar have noticed that the arithmetic branching program for 
the determinant can be transformed directly into a symmetric determinant of dimensions 
0{n'^) with techniques similar to the ones used in this paper. We give a detailed proof in 
Subsection 3.2. 

We add that the assymptotically smallest known division-free algebraic circuits for the 
nxn determinant polynomial have size 0{n?'™) [Kaltofen 1992; Kaltofen and Villard 2004]. 
The circuits actually can compute the characteristic polynomial and the adjoint and are 
based on algebraic rather than combinatorial techniques. Weakly skew circuits of such size 
appear not to be known. 

Organization. Section 1.2 is devoted to an introduction to the algebraic complexity 
theoretic used in our constructions, as well as a reminder of the existing related constructions 
in algebraic complexity. Section 2 deals with symmetric representations of formulas while 
Section 3 focuses on weakly-skew circuits. Table 2 page 39 gives an overview of all the 
different constructions used in this paper. Section 4 then proceeds to the comparisons 
between the results obtained so far and Quarez's [2008]. The special case of fields of 
characteristic 2 is studied in Section 5. 



This bound can be found on p. 11 of their paper. 



A shorter version of this paper [Grenet et al. 2011a] has been pubhshed in Proceedings of 
STAGS 2011. It contains material from Section 3 and Section 5. 
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1.2 Known results and definitions 

In his seminal paper Valiant [1979] expressed the polynomial computed by an arithmetic 
formula as the determinant of a matrix whose entries are constants or variables. If we define 
the skinny size e of the formula as its number of arithmetic operations then the dimensions 
of the matrix are at most e + 2. The proof uses a weighted digraph construction where 
the formula is encoded into paths from a source vertex to a target, sometimes known as 
an Algebraic or Arithmetic Branching Program [Nisan 1991; Beimel and Gal 1999]. This 
theorem shows that every polynomial with a sub-exponential size formula can be expressed 
as a determinant with sub-exponential dimensions, enhancing the prominence of linear 
algebra. A slight variation of the theorem is also used to prove the universality of the 
permanent for formulas which is one of the steps in the proof of its VNP-completeness. In 
a tutorial, von zur Gathen [1987] gives another way to express a formula as a determinant: 
his proof does not use digraphs and his bound is 2e -|- 2. Refining his techniques, Liu and 
Regan [2006] gave a construction leading to an upper bound of e -|- 1 in a slightly more 
powerful model: multiplications by constant are free and do not count into the size of the 
formula. 

Our purpose here is to express a formula as a determinant of a symmetric matrix. 
Multiplications by constant are also given for free. Our construction uses paths in graphs, 
similar to the paths in digraphs in Valiant's original proof. In fact, this original construction 
appears to have a little flaw in it. Interestingly enough, this flaw has never been mentioned 
in the literature to the authors' knowledge. A slight change in the proof is given in [Biirgisser 
et al. 1997, Exercise 21.7 (p570)] that settles a part of the problem. And the same flaw 
appears in the proof of the universality of the permanent in [Biirgisser 2000] . When adding 
two formulas, the resulting digraph can have two arcs between the source and the target, 
which can lead to the sum of two variables being an entry of the matrix, and this is not 
allowed as we seek for symmetric matrices where each entry is either a constant or a variable. 
The first idea to correct the proof is to keep the same parity for all s-t-paths as in Valiant's 
original proof, adding two new vertices and replacing one of the arcs by a length-three path. 
This method is very simple but its disadvantage is that it increases the dimensions of the 
final matrix to 2e -|- 3. In the symmetric case we will use a —1 coefficient to correct the 
parity differences between paths instead of adding new vertices. Using this technique in 



the non-symmetric case allows us to prove Valiant's theorem with (e + 1) instead of (e + 2). 
Our technique also gives for free multiplications by constants as in [Liu and Regan 2006]. 
It uses digraphs and is to our opinion more intuitive than direct work on matrices. 

In [Toda 1992; Malod and Portier 2008], results of the same flavor were proved for a 
more general class of circuits, namely the weakly-skew circuits. Malod and Portier [2008] can 
deduce from those results a fairly simple proof of the VQP-completeness of the determinant 
(under g^p-projection). Moreover, they define a new class \/Pws of polynomials represented 
by polynomial-size weakly-skew circuits (with no explicit restriction on the degree of the 
polynomials) for which the determinant is complete under p-projection. A formula is a 
circuit in which every vertex has out-degree 1 (but the output). This means in particular 
that the underlying digraph is a tree. A weakly-skew circuit is a kind of generalization of 
a formula, with a less constrained structure on the underlying digraph. For an arithmetic 
circuit, the only restriction on the digraph is the absence of directed cycles (that is the 
underlying digraph is a directed acyclic graph). A circuit is said weakly-skew if every 
multiplication gate a has the following property: the sub-circuit associated with one of its 
arguments /3 is connected to the rest of the circuit only by the arrow going from (3 to a. 
This means that the underlying digraph is disconnected as soon as the multiplication gate 
a is removed. In a sense, one of the arguments of the multiplication gate was specifically 
computed for this gate. 

Toda [1992] proved that the polynomial computed by a weakly-skew circuit of skinny 
size e can be represented by the determinant of a matrix of dimensions (2e -|- 2). This 
result was improved by Malod and Portier [2008]: The construction leads to a matrix of 
dimensions (m -|- 1) where m is the fat size of the circuit (i.e. its total number of gates, 
including the input gates). Note that for a circuit in general and for a weakly-skew circuit 
in particular m < 2e + 1. The latter construction uses negated variables in the matrix. It 
is actually possible to get rid of them [Kaltofen and Koiran 2008[. Although the skinny 
size is well suited for the formulas, the fat size appears more appropriate for weakly-skew 
circuits. In Section 3, we symmetrize this construction so that a polynomial expressed by a 
weakly-skew circuit equals the determinant of a symmetric matrix. Our construction yields 
a symmetric matrix of dimensions (2m + 1). In fact, this can be refined as well as the 
non-symmetric construction. An even more appropriate size for a weakly-skew circuit is 
(e -|- i) where e is the skinny size and i the number of inputs labelled by a variable (clearly 
e + i < m). We can show that the bounds are still valid if we replace m by (e -|- i) and even 
when multiplications by constants are free as in [Liu and Regan 2006] (see Section 3.3). 

Let us now give some formal definitions of the arithmetic circuits and related notions. 

Definition 1. An arithmetic circuit is a directed acyclic graph with vertices of in-degree 
or 2 and exactly one vertex of out-degree 0. Vertices of in-degree are called inputs and 
labelled by a constant or a variable. The other vertices, of in-degree 2, are labeled by x 
or -f and called computation gates. The vertex of out-degree is called the output. The 
vertices of a circuit are commonly called gates and its arcs arrows. 



An arithmetic circuit with constant inputs in a field k and variables in a set x naturally 
computes a polynomial / G k[x]. 

Definition 2. If q is a gate of a circuit C, the sub-circuit associated to a is the subgraph of 
C made of all the gates /3 such that there exists a oriented path from /3 to a in C, including 
a. A gate a receiving arrows from /3 and 7 is said to be disjoint if the sub-circuits associated 
to P and 7 are disjoint from one another. The gates /? and 7 are called the arguments of a. 

Definition 3. An arithmetic circuit is said weakly-skew if for any multiplication gate a, 
the sub-circuit associated to one of its arguments /3 is only connected to the rest of the 
circuit by the arrow going from /3 to a: it is called the closed sub-circuit of a. A gate which 
does not belong to a closed sub-circuit of C is said to be reusable in C. 
A formula is an arithmetic circuit in which all the gates are disjoint. 

The reusability of a gate depends of course on the considered circuit C. For instance, 
in Fig. 1(b), the weakly-skew circuit has two closed sub-circuits. The input z is in the 
right closed sub-circuit and is therefore not reusable. But inside this closed sub-circuit, it 
is reusable, and actually used as argument to the summation gate twice. Figures 1(a) and 
(c) are respectively an equivalent arithmetic circuit and an equivalent formula, that is the 
two circuits and the formula compute the polynomial {x -\- y)^ -|- 2yz. 

Let us remark a fact that will be useful later: all the multiplication gates of a weakly- 
skew circuit are disjoint (but this is not a sufficient condition). 

In our constructions, we shall use graphs and digraphs. In particular, the improved 
construction based on Valiant's represents formulas by paths in a digraph. On the other 
hand, to obtain symmetric determinantal representations the digraphs have to be symmet- 
ric. These correspond to graphs. In order to avoid any confusion between directed and 
undirected graphs, we shall exclusively use the term graph for undirected ones, and other- 
wise use the term digraph. It is well-known that cycle covers in digraphs are in one-to-one 
correspondence with permutations of the vertices and therefore that the permanent of the 
adjacency matrix of a digraph can be defined in terms of cycle covers of the digraph. Let 
us now give some definitions for those facts, and see how it can be extended to graphs. 

Definition 4. A cycle cover of a digraph G = {V, A) is a set of cycles such that each vertex 
appears in exactly one cycle. The weight of a cycle cover is defined to be the product of 
the weights of the arcs used in the cover. Let the sign of a vertex cover be the sign of the 
corresponding permutation of the vertices, that is (—1) where A^ is the number of even 
cycles. Finally, let the signed weight of a cycle cover be the product of its weight and sign. 
For a graph G = (V, E), let G = {V, A) be the corresponding symmetric digraph. Then 
a cycle cover of G is a cycle cover of G'^, and the definitions of weight and sign are extended 
to this case. In particular, if there is a cycle cover of G with a cycle G = (ui, . . . ,Uk), 
then a new cycle cover is defined if G is replaced by the cycle (n^, . . . ,ui). Those two cycle 
covers are considered as different cycle covers of G. 




(a) An arithmetic 
circuit 




(b) A weakly-skew circuit: The 
closed sub-circuits are marked with 
gray boxes. 




(c) A formula 

Figure 1: An arithmetic circuit, a weakly-skew circuit and a formula computing the same 
polynomial (x + y)'^ + "^yz. 




Figure 2: A graph with a cycle cover (the arcs in the cover are in black). 



Definition 5. Let G be a digraph. Its adjacency matrix is the (n x n) matrix A such that 
Aij is equal to the weight of the arc from i to j {Aij = is there is no such arc). The 
definition is extended to the case of graphs, seen as symmetric digraphs. In particular, the 
adjacency matrix of a graph is symmetric. 

Lemma 1. Let G be a (di) graph, and A its adjacency matrix. Then the permanent of A 
equals the sum of the weights of all the cycle covers of G, and the determinant of A is equal 
to the sum of the signed weights of all the cycle covers of G. 

Proof. The cycle covers are obviously in one-to-one correspondence with the permutations 
of the set of vertices, and the sign of a cycle cover is defined to match the sign of the 
corresponding permutation. Suppose that the vertices of V are {1, . . . ,n} and let Aij be 
the weight of the arc (i,j) in G. Let C a cycle cover and a the corresponding permutation. 
Then it is clear that the weight of G is Ai^^ii\ ■ ■ ■ A^^^i^\, hence the result. D 

The validity of this proof for graphs follows from the definition of the cycle covers of a 
graph in terms of the cycle covers of the corresponding symmetric digraph. In the following, 
the notion of perfect matching is used. A perfect matching in a graph G is a set M of edges 
of G such that every vertex is incident to exactly one edge of M. The weight of a perfect 
matching is defined in this as the weight of the corresponding cycle cover (with length-2 
cycles). This means that this is the product of the weights of the arcs it uses, or equivalently 
it is the square of the product of the weights of the edges it uses. Note that this is the 
square of the usual definition. 

A path P in a digraph is a subset of vertices {ui, . . . , Uk} such that for 1 < i < k — 1, 
there exists an arc from Ui to Uj+i with nonzero weight. The size \P\ of such a path is k. 

2 Formulas 

2.1 Non-symmetric case 

In this section, as in Sections 2.2 and 3, a field k of characteristic different from 2 is fixed 
and the constant inputs of the formulas and the weakly-skew circuits are taken from k. The 
variables are supposed to belong to a countable set x = {xi,X2, . . . }. Following [Liu and 
Regan 2006], we define a formula size that does not take into account multiplications by 
constants. 

Definition 6. Consider formulas with inputs being variables or constants from k. The 
green size gsize((/?) of a formula 93 is defined inductively as follows: 

• The green size of a constant or a variable is 0; 

• If c is a constant then the green size oi ex ip \s equal to the green size of <y9; 

• If (fi and 992 are formulas, then gsize((/?i -|- (^2) = gsize((/?i) -|- gsize((/?2) + 1- 



• If 991 and ip2 are non-constant formulas, then gsize(c/9i x (^2) = gsize(c/9i)+gsize((^2) + l 

An even smaller size can be defined by deciding that every variable-free formula has size 
zero and Theorem 1 can easily be extended to this case. A formal definition of this size is 
given is Section 3.3 in the context of weakly-skew circuits. 

Theorem 1 ([Liu and Regan 2006]). For every formula 99 of green size e with at least one 
addition there is a square matrix A of dimensions e + 1 whose entries are inputs of the 
formula and elements of {0, 1, —1, 1/2} such that f = det{A). 

We remark that if ip has no addition it is of the form cxi . . . x„ and it has size (n — 1). 
Then a suitable matrix is the (n -|- 1) x (n + 1) diagonal matrix made of the n variables 
and the constant c. Thus the dimensions of the matrix are at most n + 1 = e + 2, and are 
?i = e + lifc=l. Note that this latter bound is minimal as the determinant of a (d x d) 
matrix is a degree-d polynomial. The dimensions (n + 1) are not minimal when c 7^ 1 as 
shown by the (3 x 3) matrix 
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representing 2xyz. One can also see that the n bound cannot be general as there is no 
(2 X 2) matrix representing the polynomial 2xy. 

Lemma 2. Let ip be an arithmetic formula of green size e. Then there exists a constant cq 
and an edge-weighted digraph G with at most e + 2 vertices and two distinct vertices s and 
t such that 

s-t-path P 

Proof of Lemma 2. We prove the lemma by induction on formulas. If 93 is equal to a variable 
X (resp. a constant c) then G has two vertices s and t and an edge (s, t) labelled by x (resp. 
c) and the constant cq is equal to 1. 

If 99 = c X 99' let G' be the digraph and Cq the constant satisfying the lemma for the 
formula ip' . Then obviously G = G' and cq = CgC satisfy the lemma for 93. 

\i Lp = ipi X (p2, let Gi and ci (resp. G2 and C2) satisfying the lemma for ipi (resp. 992). 
Then let c = C1C2 and G be the disjoint union of Gi and G2, except for ti and S2 which are 
merged (see Fig 3). The size of G is equal to |Gi | -|- IG2I — 1 < gsize(</?i) -|- gsize(992) + 3 = 
gsize(99) + 2. A s-t-path P in G is a si-ti-path Pi in Gi followed by a S2-t2-psith P2 in G2 
and we have |P| = |Pi| -|- IP2I — 1 and w{P) = w{Pi) x w{P2), hence the result. 

li (p = ipi -\- ip2, let Gi and ci (resp. G2 and C2) satisfying the lemma for 991 (resp. 992). 
If ci =0 then ip and 992 compute the same polynomial and we just have to take G = G2 and 
c = C2- Suppose now ci 7^ 0. Then we define G as the disjoint union of Gi and G2, except 
for si and S2 which are merged, and with an edge (^25^1) of weight — C2/C1 (see Fig 4). The 
size of G satisfies the same relation as in the multiplication case. Let cq = Ci. A s-t-path P 
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Figure 3: Gi,ci and G2,C2 are respectively associated to (pi and ip2; f = (fi x (p2- 
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Figure 4: Gi,ci and G2,C2 are respectively associated to c^i and 992; (/? = c^i + (/72 • 



in G is a si-ti-path in Gi or a S2-t2-path P2 in G2 followed by the edge {t2,ti), and in the 
second case we have ■w{P) = w{P2){—C2/ci) and \P\ = IP2I + !> hence the result. Remark 
that ^2 has only one outgoing edge and its weight is a constant, and that this property will 
not be changed in the inductive construction. This property will be useful to prove the 
bound in the theorem. D 



Proof of Theorem 1. Let ^p be an arithmetic formula of green size e and let G and cq be 
given by Lemma 2. Let G be the digraph obtained from G in the following way. We merge 
s and t. As remarked in the proof of Lemma 2 there is a vertex v that has only one outgoing 
edge and its weight is a constant c (as (^ is supposed to have at least one addition). We 
change its weight to cqc and add a loop weighted by cq on v. We put a loop with weight 1 
on every other vertex than v and s. 

Let {1, . . . , e + 1} be the vertices of G and A its adjacency matrix. Let us have a closer 
look at cycle covers of G. The cycles in G are cycles containing s (which are in bijection with 
■s-t-paths in G) and loops. In a cycle cover G the vertex s belongs to a cycle S. Its weight 
w{s) is the weight of the corresponding s-t-path P in G and its cardinal is 15*1 = \P\ — 1. 
If the vertex v appears in S then w{S) = cqw{P) and every other cycle in G is a loop of 
weight 1. Otherwise w{S) = w{P) and G contains the loop v of weight cq. In both case 
w{G) = cqw{P). Let us recall that sgn(G) is the signature of the underlying permutation: 
here it is —1 if S is even and 1 otherwise, and so it is equal to (— 1)' L Using Lemma 1 we 
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get 

det(^)= Yl sgn(CMC) = co- Yl ("1)"'' ^(^) = <^- 

cycle cover s-t-path 

C of G ^6G 



D 



2.2 Symmetric case 

The aim of this section is to write an arithmetic formula as a determinant of a symmetric 
matrix, whose entries are constants or variables. Recall that in this section as in Section 3, 
a field k of characteristic different from 2 is fixed, and the input constants are taken from 
this field. In the sequel, every constructed graph is undirected. At first, the result is proved 
for the skinny size of the formula. We recall that the skinny size of y? is the number of 
arithmetic operators it contains. 

Theorem 2. Let (p he an arithmetic formula of skinny size e. Then there exists a matrix 
A of dimensions at m,ost 2e + 3 whose entries are inputs of the formula and elements of 
{0, 1, -1, 1/2} such that ip = det A. 

This theorem is a corollary of the following lemma. 

Lemma 3. Let ip he an arithmetic formula of skinny size e. Then there exists a graph G 
with at most 2e + 2 vertices and two distinct vertices s and t such that 

1. The graph G has an even number of vertices, every cycle in G is even and every 
s-t-path has an even number of vertices. 

2. The subgraph G\ {s,t} is empty if e = and for e > 1 it has only one cycle cover: It 
is a perfect matching of weight 1. For every s-t-path P in G, the subgraph G\P is 
empty or has only one cycle cover: as above it is a perfect matching of weight 1. 

3. The following equality holds in G: 



Y r-nl^l/2+1 



s-t-path P 

The graph G is called the graph associated to ip. 

The first property of the lemma ensures that because of a parity argument every cycle 
cover of the final constructed graph G used in the proof of Theorem 2 (see Fig. 5) includes 
exactly one path between s and t. The second property ensures that the weight of the cycle 
cover is the weight of the cycle involving s and t, that is every other cycle has weight 1, 
and that other cycles of the cover are of length 2. The third property gives the relation 
between the graph and the formula. 

11 



As in Valiant's construction for the not necessarily symmetric case, the formula 99 will 
be encoded in the weights of paths between s and t, but in a slightly different way. In 
Valiant's construction, a cycle cover of the digraph is made of a cycle including a s-i-path, 
other cycles being loops. Moreover every s-t-path has the same parity and so every cycle 
cover has the same parity of odd cycles and the underlying permutation has the same 
signature. With this property of the digraph the determinant of its adjacency matrix is 
equal to its permanent up to the sign. In our construction a cycle cover of the graph is 
made of a cycle including a s-i-path, other cycles being length-2 cycles. A length-2 cycle 
has a negative signature and every s-i-path of the graph has an even cardinality, so the 
sign of the cycle permutation is —1 to the number of length 2 cycles. This shows that the 
sign of the cycle permutation is a function of the length of the involved s-t-path modulo 4. 
There is a way to ensure that this sign does not depend on the chosen s-t-path: replace the 
graph G associated to a size-0 formula x in the proof of Lemma 3 by a 4-vertices path with 
weight X on its first edge, and replace weights —1 (Fig. 5, Fig. 7 and Fig. 8) by weights 1. 
This yields a matrix with entries in kU x whose determinant and permanent are equal to 
(p, but its dimensions can be 4e + 5. To achieve the 2e + 3 bound, we construct a matrix A 
whose determinant can be very different from the permanent: For example, the permanent 
of the matrix associated to 93 = x + 2; is when its determinant is 2x. Nonetheless we can 
very easily obtain a matrix B having the same dimensions as A and such that perm B = ip 
by replacing every —1 entry in A by 1. 



Proof of Theorem 2. Let G be the graph associated to p and let G be the graph G aug- 
mented with a new vertex c and the edges tc of weight 1/2 and cs of weight ( — l)l'^l/2-i 
(see Fig. 5). 




Figure 5: Construction of G from G. 

Conditions (1) and (2) imply that there is a bijection between paths from s to t or t 
to s and cycle covers in G. More precisely, every cycle cover in G has a unique odd cycle 
and it is of the form cPc where P is a s-t-path or a t-s-path. Indeed, the graph G has an 
odd number of vertices. Suppose there is a cycle cover of G involving the length-2 cycle 
tct. Other cycles of this cover are cycles of G and thus by (1) they are all even. This is 
not possible as an odd set can not be partitioned into even subsets. For the same reason, 
there is no cycle cover of G involving the cycle scs. Thus every cycle cover of G has a cycle 
including c and a path P between s and t. 
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Let us recall that the sign of a cycle cover is the sign of the underlying permutation, 
i.e. —1 if it has an odd number of even cycles and 1 otherwise, and let us define the signed 
weight of a cycle cover as the product of its weight and sign. Let C be a cycle cover of G 
involving the s-i-path P. By property (2) there is only one way to complete the cover. Thus 
the weight of the cycle cover is the weight of P multiplied by (1/2 (— l)l'^l/2+i) and its sign 
is the sign of a perfect matching of cardinality |G \ P|, so it is ( — l)(l'^\-P|)/2, By symmetry, 
the inverse cycle cover has the same signed weight. So the sum of the signed weights of all 
cycle covers of G is equal to twice the sum over all s-t-path P of (1/2 ( — 1)1^1/2+1 w{P)). 
According to Lemma 3 it is equal to y?. The result follows from Lemma 1. 

D 

Proof of Lemma 3. We proceed by structural induction. In other words, we first prove the 
lemma for the simplest possible formula, namely x, and then show that the assertion of 
Lemma 3 is stable under addition and multiplication. 

Let (/9 = X be an arithmetic formula of size 0. Then the graph G associated to ip by 
definition has two vertices s and t and an edge st of weight x. It verifies trivially properties 
(1) and (2) and its only s-t-path is st and we have: ( — l)^'^"'""^x = ip. 

Let if = ipi + ip2 and Gi and G2 be the graphs associated to ipi and c/?2- First let us 
suppose siti or ^2^2 has weight 0. This means in particular that (pi or ip2 is of size at least 1. 
Let s = si = S2 and t = ti = ^2- Suppose Gi\{si,ti} and G'2\{s2,t2} have disjoints sets of 
vertices and let G = GiUG2 (see Fig. 6). Then |G| = 1^11 + 1021-2 < 2\ipi\+2\ip2\+2 = 2|(^|. 






^2 





© 



Figure 6: Graph associated to ip = ipi + ip2. 



If siti is an edge in Gi and S2t2 is an edge in G2 then the preceding construction 
would lead to two edges between s and t. They could be transformed into a single edge 
if adding the two weights, but then the weight could be a sum of two variables, and it is 
something that is not allowed in this context. So the graph Gi is transformed into a graph 
G'l by adding two vertices u and v, removing the edge siti with weight x and adding the 
edges siu with weight x, uv with weight 1 and vti with weight —1 (see Fig. 7). We can 
verify easily that G'l satisfies the three conditions of Lemma 3. In particular for the third 
condition, the term x corresponding to the path siti in Gi in the sum is replaced by the 
term corresponding to the path s\uvti in G'^: —{—l)^''^~^^x = x. We then construct the 
graph G associated to ip as above but with G'l replacing Gi. It size is at most 2\ip\ + 2. 
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Figure 7: Transformation of Gi into G'^. 



Now let us prove that the graph associated to (p satisfies the three properties of the 
lemma. 

1. G has an even number of vertices and the cardinality of every s-t-path is even. A 
cycle in G is a cycle in Gi , or a cycle in G2 , or a path from s to t in Gi or G2 followed 
by path from i to s in Gi or G2, and consequently every cycle in G is even. 

2. If Gi \ {si,ti} and G2 \ {s2,t2} are non-empty they are disconnected, and a cycle 
cover of the subgraph G \ {s, t} is constituted by a cycle cover of Gi \ {si, ti} and a 
cycle cover of G2 \ {525^2}- So G \ {s,t} has only one cycle cover and it is a perfect 
matching of weight 1. If Gi \ {si,ti} is empty then G \ {s, t} = G2 \ {521^2} and has 
only one cycle cover and it is a perfect matching of weight 1 . 

Let P be a path between s and t in G. We can suppose wlog that the subgraph 
G \ P is the union of the two graphs Gi\P and G2 \ {s2, ^2}, which are disconnected 
from one another. The property to prove is then straightforward from the induction 
hypothesis. 

3. A path of G is a path of Gi or a path of G2, which proves the equality. 

Let ip = ipiXip2 and Gi and G2 be the graphs associated to ipi and ip2- Suppose Gi and 
G2 have disjoints sets of vertices and let G be Gi UG2 with an additional edge tiS2 of weight 
-1, and let s = si and t = t2 (see Fig.8). Then \G\ = |Gi[ + IG2I < 2\ipi\ + 2\ip2\ + 4 = 
2\ip\ + 2. Let us prove that G satisfies the three properties of the lemma. 

1. G has an even number of vertices and every path from s to t has an even cardinality. 
A cycle in G is either a cycle in Gi, or a cycle in G2 or the length-2 cycle tiS2, and 
consequently every cycle in G is even. 

2. Let us consider a cycle cover of G \ {s,i}. The vertex ti can be in a cycle of Gi or 
in the cycle tiS2- If it is in a cycle of Gi then we have a cycle cover of Gi \ {si}, 
which is not possible because it is an odd set and all its cycles are even. Thus the 
cycle cover of G \ {s, t} can be partitioned into tiS2 of weight ( — 1)^, a cycle cover of 
Gi \ {si,ii} and a cycle cover of G2 \ {52,^2}- Those cycle covers are unique and so 
there is only one cycle cover of G \ {x, y} and it is a perfect matching of weight 1. 
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Figure 8: Graph associated to 99 = 991 x (p2. 



Let P be a path between s and t in G. It is a path Pi from si to ti in Gi followed 
by tiS2 and a path P2 from S2 to ^2 in G2. So G \ P is the union of the two graphs 
Gi \ Pi and G2 \P2, which are disconnected (if non empty) from one another. The 
property to prove is then straightforward from the induction liypothesis. 

3. A s-t-path P in G can be decomposed into three paths: a si-ti-path Pi, tiS2 which 
is of weight — 1 and a S2-i2-path P2 . 

Thus 



(-1) 



\p\ 



+1 



w{P) 



(_1)^+1 w{Pi) X (-1)^+1 W{P2) 



and so 



E 



\p\ , 1 
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Pi 

<fl X (p2 
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-P2 
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D 

The upper bound (2e + 2) of Lemma 3 is tight as shown by Fig. 9. It can be shown 
easily that this construction yields a graph of size at least \ip\ +2, and this lower bound is 
tight as shown by Fig. 10. 

In fact, as in the non-symmetric case, the skinny size can be replaced by the green size 
of the formula defined in Definition 6. 

Theorem 3. For every formula ip of green size e there is a square matrix A of dimensions 
2e + 3 whose entries are inputs of the formula and elements of {0, 1, —1, 1/2} such that 
if = det A. 
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Figure 9: Graph associated to ip = xi + ■ ■ ■ -\- Xn+i- \^p\ = n and |G| = 2n + 2. 




Figure 10: Graph associated to ip 



Xix'i +X2X2 ' 



■ + XnXn + y: \ip\ 



2nand \G\ = 2n + 2. 



Proof. It is sufficient to show how to have the constants for free in the construction of 
Lemma 3. We also proceed by structural induction. In fact, the construction remains 
almost the same but with the last property changed. For an arithmetic formula if of green 
size e, there exists a graph G that satisfies the conditions of Lemma 3 but the third one is 
replaced by the existence of a constant cq such that 

CO- Y. i-^r^'^' HP) = ^. 

s-t-path P 

Let c^ = X be an arithmetic formula of size 0. Then the graph G associated to tp by 
definition has two vertices s and t and an edge st of weight x. The associated constant is 
Co = 1. 

Let (f = cip and G, cq be associated to ip. Then G, ccq is associated to ip. 

Let ip = (pi X ip2 and Gi, ci (resp. G2, C2) be associated to ipi (resp. 1^2)- The graph 
G associated to ip is exactly the same as in the proof of Lemma 3 and the constant is C1C2. 

Let ip = ipi + {p2 and Gi, ci (resp. G2, C2) be the graph and constant associated to ipi 
(resp. ip2). We suppose that Gi and G2 have distinct sets of vertices except for si = 32- 
The graph G is obtained by adding a new vertex u, an edge t2U with weight 1 and an edge 
uti with weight — C2/C1, and the associated constant is ci (see Fig. 11). 

This defines a size-(2e + 2) graph G associated to a green size-e formula ip. It remains 
to turn this graph into a matrix. Let G be the graph G augmented with a new vertex c 
and the edges tc of weight Co/2 and cs of weight {—ly'-''''^~^. The adjacency matrix A of 
G satisfies 99 = det(^) and the proof is similar to the one of Theorem 2. D 
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Figure 11: ip = (pi + ip2', Gi,ci and G2, C2 are respectively associated to ipi and ip2- 



The bound obtained in Theorem 3 can be sharpened when /c = M or C The idea is to 
build G by merging s and t instead of adding a new vertex. Suppose that 93 has at least one 
addition gate. Let w = y^\co\/2. In the construction for this addition gate (see Fig. 11), 
multiply the weights of t2U and uti by w. A cycle cover of the graph either goes through 
the path t2uti, or contains the edge ut2 in its perfect matching part. In both cases, its 
weight is multiplied by w"^. Now if {—iy'^'''^~^^co/2 > 0, then the graph obtained has the 
satisfying properties, and the new bound is 2e + 1. If it is negative, two solutions can be 
applied. Either k is the field of complex numbers and it is sufficient to replace w by iw 
(where i'^ = —1) to get the same bound 2e + 1. Otherwise, if k is the field of real numbers, 
it is sufficient to add a new vertex with a loop of weight —1 (this corresponds to adding a 
new line and a new column, filled with zeroes but the diagonal element with —1) to get the 
bound (2e + 2). 

3 Weakly skew circuits 

In this section, we extend the previous results to the case of weakly-skew circuits. Recall 
that those circuits are defined from arithmetic circuits by a restriction on the multiplication 
gate: the sub-circuit associated to one of the arguments of a multiplication gate a has to 
be closed, that is only connected to the rest of the circuit by the arrow going to a. A gate 
that is not in any such closed sub-circuit is said to be reusable. 

The main difficulty to extend the results is the existence of several reusable gates. In 
the case of formulas, there is a single output. Therefore, there is a single vertex t in the 
graph for which the sum of the weights of the s-t-paths has to equal a given expression. 
This is no longer the case for weakly-skew circuits. If the matrix we wish to construct is 
not symmetric, that is if the graph is oriented, this difficulty is overcome by ensuring that 
the graph is a directed acyclic graph. In that way, adding a new vertex cannot change the 
expressions computed at previously added vertices. But in the symmetric case, adding a 
new vertex, for example in the case of an addition gate, creates some new paths in the 
graph. Thus it changes the sum of the weights of the s-t^-paths for some vertex ta- 

A solution to this problem is given in Lemma 4 by introducing the notion of acceptable 
paths: A path P in a graph G is said acceptable ii G\P admits a cycle cover. 
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3.1 Symmetric determinantal representation 

For the weakly-skew circuits, the green size is no longer appropriate. Hence, the results of 
this section are expressed in terms of the fat size of the circuits: the fat size of a circuit is 
its total number of gates, including the input gates. This measure of circuit size is refined 
in Section 3.3. 

Theorem 4. Let f be a polynomial computable by a weakly-skew circuit of fat size m. Then 
there exists a symmetric matrix A of dimensions at most 2m + 1 whose entries are inputs 
of the circuit and elements from {0, 1, —1, 1/2} such that f = det j4. 

The proof relies on the following lemma. It applies to so-called multiple- output weakly- 
skew circuits. This generalization just consists of circuits for which there exist several 
out-degree-0 gates. 

Lemma 4. Let C be a multiple- output weakly- skew circuit of fat size m. There exists a 
graph G with at most 2m -\- 1 vertices and a distinguished vertex s such that \G\ is odd, 
every cycle in G is even, and for every reusable gate a G C there exists a vertex to, G G 
such that 

1. Every s-ta-path (whether acceptable or not) has an odd number of vertices; 

2. For every acceptable s-ta-path P in G, the subgraph G \ P is either empty or has a 
unique cycle cover, which is a perfect matching of weight 1; 

3. The following equality holds in G: 

^ (-l)^^P) = /„ (2) 

acceptable 
s-ta-path P 

where f^ is the polynomial computed by the gate a. 

Furthermore, the graph G \ {s} has a unique cycle cover which is a perfect matching of 
weight 1. 

Proof. The graph G is built by induction on the (fat) size of the circuit, the required 
properties being verified at each step of the induction. If a is a reusable gate of C, then ta 
is said to be a reusable vertex of G. 

A size-1 circuit is an input gate a with label x. The corresponding graph G has three 
vertices: s, ta and an additional vertex Va- There is an edge between s and Va of weight x, 
and an edge between Va and ta of weight —1. It is straightforward to check that G satisfy 
the conditions of the lemma. 

Let m > 1 and suppose that the lemma holds for any multiple-output weakly-skew 
circuit of size less than m. Let C be a multiple output weakly-skew circuit of size m, and 
a be any of its outputs. 
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If a is an input gate with label x, let C' = C \ {a} and G' the corresponding graph 
with a distinguished vertex s. The graph G is obtained from G' by adding two new vertices 
Va and ta, an edge of weight x between s and Va and an edge of weight —1 between Va 
and ta (see Fig. 12). The vertex s is the distinguished vertex of G. The size of G is 



c 



Figure 12: Induction step when a is an input gate. 

|G| = \G'\ + 2 < (2(m - 1) + 1) + 2 = 2m + 1. Thus |G| is odd. A cycle in G is either a 
cycle in G' or one of the two cycles sVa or Vata, so every cycle in G is even. The size-3 path 
from s to ta is acceptable (as G' \ {s} has a unique cycle cover of weight 1) and satisfies (2). 
Now, any other reusable gate /3 belongs to G', so the conditions are satisfied by induction 
hypothesis (it is sufficient to remark that when s is removed, Va and ta are disconnected 
from the rest of the circuit, and a cycle cover has to match those two vertices). 

If a is an addition gate, let G' = G \ {a} and suppose that a receives arrows from 
gates /3 and 7. Note that (3 and 7 are reusable. Let G' be the graph corresponding to C, 
and s be its distinguished vertex. G' contains two reusable vertices i^ and t^. The graph 
G is obtained by adding two vertices Va and ta, and the following edges: tj^Va and t^Va 
of weight 1, and Vata of weight —1 (see Fig. 13). If /3 = 7, then G' contains a vertex i^, 
and we merge the two edges adjacent to t^ and t^ into an edge tpVa of weight 2. Then 



; G' 








Figure 13: Induction step when a is an addition gate. 



|G| = |G'| + 2 < 2m + 1, and \G\ remains odd. 

Every s-t^-path for some reusable gate 5 in G' is even. A cycle in G is either a cycle 
in G' , or the cycle Vata, or is made of a t^-t^-path P in G' plus the vertex Va- Let P' 
be a s-t/3-path and u the first vertex of P' belonging to P. Then, P' = s, . . . ,u, . . . ,t/3 
and P" = s, . . . ,u, . . . ,t:y are both path with an odd number of vertices. In particular the 
sizes of u, . . . , t^ and u, . . . ,t^ are of same parity. Thus P is of odd size and P U {va} is 
an even-size cycle. Hence, every cycle in G is even. An acceptable path in G is either an 
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acceptable path in G' or a path from s to ta- Indeed, the only way to cover ta in a cycle 
cover is to match it with v^. Therefore, no acceptable path goes through tp, Va and t-y. So, 
the reusable gates in C' satisfy the conditions of the lemma by induction. Any acceptable 
path P from s to ta is an acceptable path P' from s to t^ or ty followed by a path from t^ 
or t^ to ta- Thus \P\ = \P'\ + 2 is odd and G\P = G' \P' has a unique cycle cover which 
is a perfect matching of weight 1. Finally, 



acceptable 
s-ta-path P 



|f/3l+2-l 



E i-l)^^ i-l ■ w{Pp)) + E (-l)^^^(-l-^^7)) 



acceptable 
s-t^-path P0 



acceptable 
s-i^-path Pry 



|f/3l-l 



Pj\-1 



^i-l)^^ w{Pfs) + E(-l)^^^7) 



—fl3 + f-y — fa- 

If a is a multiplication gate, a receives arrows from two distinct gates /3 and 7. Exactly 
one of those gates, say /3, is not reusable and removing the gate a yields two disjoint circuits 
Ci and C2 (say /3 belongs to Ci and 7 to C2). Let Gi and G2 be the respective graphs 
obtained by induction from Ci and C2, with distinguished vertices si and S2 respectively. 
The graph G is obtained as in Fig. 14 as the union of Gi and G2 where t^ and si are merged, 
the distinguished vertex s of G being the distinguished vertex S2 of G2 , and ta being equal 
to t^. Then |G| = |Gi| + IG2I — 1, so |G| is odd, and if mi and 7712 are the respective sizes 



Ci 



Co 



M?^^ '-@-'' 




; Go 1 



a — tp 



Figure 14: Induction step when a is a multiplication gate. 

of Ci and C2 {m = m\ + m2 + 1), then |G| < 2mi + 1 + 2m2 + 1 — 1 = 2m — 1. A cycle 
in G is either a cycle in Gi or a cycle in G2 and is therefore even. The reusable gates of G 
are a and the reusable gates of G2 (by definition, Gi is closed and in particular tp is not 
reusable). A path (in G) from s to a reusable gate of G2 cannot enter Gi so the reusable 
gates of G2 satisfy the first and the third conditions in the lemma. Furthermore, if such a 
path P is removed from G, the only cycle cover of G\P has to be made of a cycle cover of 
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G2\P and a cycle cover of Gi \si. Indeed, the vertex si = t^ has to be either in a cycle 
cover of Gi or in a cycle cover of G2. But G2 \ (-P U {t^}) is a graph of odd size and cannot 
be covered by cycles of even size and Gi is also of odd size. Thus, the reusable gates in G2 
also satisfy the second condition of the lemma. It remains to prove that the reusable gate 
a satisfies the conditions of the lemma: 

1. A s-ta-path P is a s-t..y-path P^ followed by a si-t^-path P/j . Thus |P| = |P.y | + |P^ | — 1 
as t^ = si and \P\ is odd. 

2. The graph G\P is the disjoint union of {G2\P'y) and (Gi \Pi3), so by induction G\P 
is either empty or has a unique cycle cover which is a perfect matching of weight 1. 

3. As w{P) = w{P-f)w{Pi3), we have 

{-l)^w{P) = (-1) -2 w{P^)w{P^) 

= {-l)^-^^w{P^) X {-l)-^w{P^), 
whence 

p Pj Pp 

= f-yX fl3 

Finally, the only way to cover G \ {s} is to cover G2 \ {82} on one hand and Gi \ {si} on 
the other hand for parity reasons as before. The weight of this cover is the product of the 
weights of the covers of Gi and G2 , that is 1 . D 

Proof of Theorem 4- Let G be a weakly-skew circuit computing the polynomial /, and G 
be the graph built from G in Lemma 4. The circuit G has a unique output, and there 
exists in G a vertex t corresponding to this output. Let G' be the graph obtained from G 

-, |G|-1 

by adding an edge between t and s of weight 2(~1) ^ ■ 

There is no cycle cover of G' containing the 2-cycle st. Indeed, \G' \ {s,t}\ is odd and 
G contains only even cycles. This means that a cycle cover of G' contains a cycle made of 
a s-t-path plus (t, s) or a t-s-path plus (s, t). Let P be such a path. Then G' \P = G\P. 
Hence, by Lemma 4, there is exactly one cycle cover of G' \P and it is a perfect matching 
of weight 1. This means that there is a one-to-one correspondence between the cycle covers 
of G' and the paths from s to t or from i to s. There is also a one-to-one correspondence 
between the paths from s to t and the paths from t to s. 

Let us recall that the sign of a cycle cover is the sign of the underlying permutation 
and its signed weight is the product of its sign and weight. Let G be a cycle cover of 
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G' involving the s-t-path P. The previous paragraph shows that the weight of C equals 

IC?I — 1 

2( — 1) 2 w{P). As C has an odd cycle and a perfect matching, its sign is (— l)l'^\^l'^, 
that is the number of couples in the perfect matching. The inverse cycle cover C of G' has 
the same signed weight as C. Hence the sum of the signed weights of all cycle covers of G' 

, \G\~1 \G\P\ -, P|-l 

equals twice the sum over all s-i-paths P of 2(~1) ^ (~1) ^ w{P) = 2(~1) ^ w{P). 
By Lemma 4, this equals / and Lemma 1 concludes the proof. 

D 

3.2 Symmetric determinantal representation of the determinant 

Let us denote by DET„ the formal determinant of the nxn matrix (xij). This polynomial 
has a weakly-skew circuit of size-O(n^) ([Berkowitz 1984; Malod and Portier 2008]) or even 
0{n'^) if we use algebraic branching program constructed by Mahajan and Vinay [1997]. 
This weakly-skew circuit can be represented by a determinant of a symmetric matrix as 
proved in this paper in Theorem 4. 

After a talk from one of us presenting our results, Meena Mahajan and Prajakta Nimb- 
horkar have communicated us the following theorem, which shows that for the determinant 
polynomial, the symmetrization can be done more efficiently that in the general case. As 
this result is not published, we find interesting to give here its proof. 

Theorem 5 (Meena Mahajan and Prajakta Nimbhorkar). For every n there is a symmetric 
matrix M of dimensions An? + 7 and entries in {xij \ 1 < i,j < n} U {0; 1; —1; 1/2} such 
that DET„ = detM. 

Proof. Construct the weighted graph G' computing DET.„ with the method used in Sec- 
tion 3 of [Mahajan and Vinay 1997]. It is a directed acyclic weighted graph with three 
distinguished vertices s, t+ and t_. Every weight is 0, 1 or a variable Xij. The graph 
satisfies 

DET„= Yl ^(^)- E ^(^) 

s—t-f-—path P in G' s—t——path P in G' 

Moreover, this graph has 2n^-|-3 vertices, at most 4n^ edges and the following nice structure: 
it is made of n-|- 1 layers, the first layer being s and the last one being {t_|_, t_}. Every edge 
is from a layer i to a layer i + 1. As a consequence, every s-t+-path has n + 1 vertices, and 
so has every s-t_-path. 

From the graph G' we can easily obtain an algebraic branching program G for computing 
DET„: add a vertex t, an edge (t+,t) of weight 1 and an edge (t-,t) of weight —1. We 
could then proceed to built from this algebraic branching program a skew-circuit of size 
0{n^) (see for example proof of Proposition 1 in [Kaltofen and Koiran 2008]) and then a 
symmetric determinantal representation of dimensions 0(?i^) with the method described 
in Theorem 4. But symmetrizing directly the algebraic branching program G allows us to 
achieve a better bound as we are going to see. 
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Let V be the set of vertices of G \ {s, t} and E be the set of edges of G. The symmetric 
weighted graph Gg is defined as fohows by duphcating vertices in graph G. The set of 
its vertices is {soutitin} U {uin,Uout \ u G ^}- The set of its edges is {uoutVin \ {u,v) G 
E} U {uinUout I u S y}- Weights on edges are defined by w{uoutVin) = w{u,v) and 
w{uinUout) = 1- The graph Gg has An^ + 6 vertices arranged in 2n + 2 layers and satisfies 
the following property: 

DET„= Y. ^(^) (3) 

acceptable 
Sout-tin-pa,th P in Ga 

Recall that a path P in a graph Gg is called acceptable if Gs\P admits a cycle cover. 

To prove Property 3, let us have a look at some acceptable Sout-^m-pa-th Pg in Gg and at 
some cycle cover C of G^ \ P^. We prove that for every u £V, the vertices iij„ and Uout are 
both in Pg or together in a length-2 cycle of C. The first vertex of the path Pg is Sout- The 
second vertex is some ni^j„ where ui is a vertex of the second layer of G. The third vertex is 
ui^out as til, in is only linked to Sgut and ui^out- Let us now consider another vertex Vin where 
V belongs to the second layer of G. It is only linked to Sout and Vout, and so it is not in Pg 
but belongs to the weight 1 and length-2 cycle VinVout in C. The same reasoning applies 
to the following layers. Thus we just proved that there is a weight-preserving bijection 
between acceptable So„i-tj„-paths in Gg and s-t-paths in G. Moreover, for every acceptable 
SoMi-^m-paths Pg in Gg, the graph Gg\P has only one cycle cover, which is of weight 1 and 
sign (— l)l<^Aftl/2 = (_i^2n -n+2 _ (—1)". Becausc of the symmetry of the graph we also 
have: 

DET„= Y. ^(^) (4) 

acceptable 
tin-Sout-pa.th P in Gs 

and thus 

DET„ = i Y HP) + l E ^(^) (5) 

acceptable acceptable 

Sout-<in-path P in Gs ti„-Sout-path P in Ga 

Remark that every cycle in Gg is even because of its layer structure. Let G be the 
graph Gg augmented with a new vertex c and the edges tj„c of weight 1/2 and cSout of 
weight (—1)", and let M be its adjacency matrix. The end of the proof is similar to the 
one of Theorem 2. The only odd cycles in G are the ones including c and a Sout-tin-psith 
or a tin-Sout-psith. P. As an odd graph can not be decomposed in even cycles, every cycle 
decomposition of in G has one of these odd cycles. It was proven above that the rest of the 
graph has only one possible cycle decomposition. Thus by (5): 

DET„= E sgn{C)w{C) (6) 

cycle cover C in G 
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According to Lemma 1 we have 

det(M) = J2 sgn(C) w{C) (7) 

cycle cover C in G 

and thus the result 

DET„ = det(M) (8) 

D 

3.3 Minimization 

The aim of this section is to refine the bound we obtained in Section 3.1, using the notion 
of green size that was defined in Section 2.1 (and matches the notion of size used in [Liu 
and Regan 2006]). As mentioned before, one can refine this notion of green size. It rehes 
on the idea ah^eady mentioned by Liu and Regan for the formulas: One can add weights on 
the arrows of the circuit. If there is an arrow from a gate a to a gate /? with weight c, then 
/3 receives as argument the value c/q where /„ is the polynomial computed by a. Such a 
circuit is called a weighted circuit. Of course, a classical circuit is a weighted one with all 
weights equal to 1. 

To refine the notion of green size, the idea is to avoid counting the variable-free sub- 
circuit. The next lemma shows that it is possible to do this in a very simple way. 

Lemma 5. If C is a weighted circuit, then there exists an equivalent weighted circuit C 
with the same number of inputs labelled by a variable and at most the same number of 
computation gates such that: 

1. An input gate is labelled either by a variable or the constant \, and the constant inputs 
have out- degree 1; 

2. An addition gate has at most one constant argument and this argument is an input 
gate; 

3. A multiplication gate has both arguments non-constant. 

Proof. One can suppose that there exists some input gate labelled by a variable, otherwise 
the polynomial computed by C would be constant. To obtain the three points, each of the 
four following rules is recursively applied to C . Each rule is applied as long as possible 
before we apply the next one. We never go back to a previous rule. 

1. Every input gate labelled by a constant c is replaced by an input gate labelled by 1, 
and the weight of an arrow going from it is multiplied by c. If there are several arrows 
going from this input gate, it is duplicated so that each copy has out-degree 1. 
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2. Every computation gate a that has both arguments constant is replaced by an input 
gate labehed by 1, and the weight of every arrow going from it is multiphed by the 
value a computed. As in previous step, the new input gates are duplicated to have 
out-degree 1. 

3. If a multiplication gate a with positive out-degree has one constant argument f3 
labelled by 1 and with an arrow from /3 to a of weight ci , and another argument 7, 
non-constant, with an arrow of weight C2, then a and /3 are deleted, and every arrow 
going from a of weight c is replaced by an arrow going from 7 of weight CC1C2 (see 
Fig. 15). 




Figure 15: Minimization for a multiplication gate. 



4. If the output gate a is a multiplication with one constant argument /3 with an arrow 
of weight ci going from /3 to a and the other argument 7, non-constant, with an arrow 
from 7 to a of weight C2, then a and /3 are deleted, 7 becomes the new output gate, 
and the weight of every arrow coming to 7 is multiplied by C1C2 (see Fig. 16). 




Figure 16: Minimization for the output gate. 



The first two rules ensure that all the constant input gates are labelled by 1 and have out- 
degree 1. After the second rule, each computation gate has at most one constant argument, 
and that it is an input gate. Then rules 3 and 4 delete all multiplication gates that have a 
constant argument. D 
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Note that the above lemma is vahd for any kind of arithmetic circuit, and that the 
construction does not change the nature of the circuit. So this can be apphed to a formula 
to get a formula, or to a weakly-skew circuit to get a weakly-skew circuit. 

Definition 7. Let C be an arithmetic circuit. Then the circuit C obtained in Lemma 5 
is the minimized circuit associated to C, and written min(C). The green size of C is equal 
to the skinny size of min(C), that is the number of computation gates in min(C). 

Note that this definition does not exactly match Definition 6 in the case of formulas, but 
is equivalent to the size mentioned right after the definition. In fact, the way of defining the 
green size we use here yields a smaller size. Nevertheless, it is easy to see that the results 
obtained in Section 2.2 remain true with this new definition. 

Theorem 6. Let f he a polynomial computable by a weighted weakly skew circuit of green 
size e and with i inputs labelled by a variable. Then there exists a symmetric matrix A 
of dimensions at most 2(e -|- i) -|- 1 whose entries are inputs of the circuit and elements of 
{0, 1, -1, 1/2} such that f = det A. 

Proof. The first step is to use Lemma 5 to minimize the circuit. Thus in the sequel the 
circuit is supposed to be a minimized weighted weakly-skew circuit. It is sufficient to show 
how to manage the constants in the construction of Lemma 4. 

The idea is to have the same construction as in Lemma 4 but with the last property 
changed: for every reusable gate a, there exists a constant Ca such that 



E (- 



■l)^-m(P) = U (9) 



acceptable 
s-ta-path P 



The changes in the construction only concern the induction steps for computation gates 
(that is for multiplication and addition gates). 

Suppose that a is an addition gate with one constant argument, say /3, with an arrow 
from /? to a of weight ci. Suppose the second argument of a is a non-constant gate 7 
with an arrow from 7 to a of weight C2. By induction, there exists a graph G^ of size 
2((e — 1) + i) + 1 that satisfies the conditions. In particular, there exists a distinguished 
vertex s, and a vertex t^ with the required properties (let c^ be the associated constant). 
Then G is obtained by adding two new vertices Va and t^ and the following edges: an edge 
t^Va of weight C2C^, an edge Vata of weight —1, and an edge sVa of weight ci (see Fig. 17). 
One can check that G satisfies the required properties. In particular, t^ satisfies (9) with 
the constant 1, and |G| = |G^| -f- 2 = 2((e - 1) + i) + 1 + 2 = 2{e + i) + 1. 

Suppose that a is an addition gate, receiving arrows from non-constant gates f3 and 7. 
There exist constants C/j and c^ such that (9) holds for f3 and 7. Suppose that the arrows 
from /3 and 7 to a have respective weights ci and C2. The construction for the induction 
step in the same as in the proof of Lemma 4, on Fig. 13, with the following changes: the 
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c\ /C2 



Figure 17: Graph obtained for the sum of a constant and a sub-circuit. 



edges tj^Va and t^Va are respectively weighted C/jCi and c^C2. Note that this does not change 
the weight of the perfect matching as those edges never belong to those matchings. As in 
that case, fa = cifp + C2f-y, we obtain 



E (-1: 



|p|-i 



acceptable 
s-ta-path P 



E {-l)^^{-ciCfs-w{Pfs)) + 



acceptable 
s-i^-path Pp 



acceptable 
s-t^-path Pj 



IP7I+2-1 



--Cl ■ Ci3 






|f;3|-l \ / |Pt|-1 \ 



=Clfl3 + C2f'y = fa- 

Note that the constant Cq associated to ta is equal to 1 in that case. If /3 = 7, with the 
same notations as above, it is sufficient to replace the weight-2 edge ti^Va by an edge of 
weight 2cic^. 

In the case of a multiplication gate, the construction (shown in Fig. 14) has no available 
edge to put the constants. But here, if the arrows from /3 and 7 to a are still labelled by 
Cl and C2 respectively, then fa = ciC2fi3f-y Thus, the same construction is kept, and the 
constant Cq associated to a is defined to be Ca = ciC2C^c-y (where c^ and c-y are respectively 
associated to /3 and 7). 

It remains to adapt the proof of Theorem 4 to this case. This is easily done by mul- 
tiplying the weight of the edge between s and t by the constant associated to the output 
gate. D 

4 Comparison with Quarez's results 

In this section, a comparison between our results and those in [Quarez 2008] is made. While 
Quarez builds matrices of fixed dimensions (depending only on the degree of the polynomial 
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and its number of variables), we build matrices whose dimensions are polynomial in the 
size of the input formula or weakly-skew circuit. Consequently, if a polynomial can be 
represented as a formula or a weakly-skew circuit of small size (say polynomial in the 
number of variables and in the degree) , then our constructions yield much smaller matrices 
than Quarez's. This is for example the case for the determinant polynomial (that is the 
determinant of a matrix of indeterminates) which is known to have a polynomial size weakly- 
skew circuit, or of the polynomial defined as the sum of all possible monomials of degree 
at most d (for this, see below). On the other hand, some polynomials are not known to 
have such polynomial size formulas or weakly-skew circuits. A famous example among 
those is the permanent. We shall see that our constructions also yield better bounds in 
that interesting case. In the most general case though, our constructions may yield bigger 
matrices. The next theorem quantifies this. 

Theorem 7. Let p be a degree-d polynomial in n variables over a field k of characteristic 
different from 2. Then p admits a formula of skinny size 

^ ' ^ - \ n+1 J \ n+1 
This yields a symmetric determinantal representation of dimensions 

5(„,d)<4("+;^-^)-2. 

Proof Let Pn^d a degree-d polynomial in n variables {xi, . . . , x„}. We shall build a weighted 
formula in the sense of Section 3.3, that is a formula with inputs in {1, xi, . . . , Xn} and with 
weights on the wires. We will first give an algorithm to build such a formula, and then 
derive an upper bound on the size of the formula so constructed. 

In order to clarify the construction, let us homogenize the polynomial Pn^d with a new 
variable xq. There exists two homogeneous polynomials Pn,d-i and Pn-i,d such that Pn,d-i 
is a polynomial of degree at most [d — 1) in (n -|- 1) variables and Pn-i,d is a polynomial of 
degree at most d in variables xq, • • • , x„_i which satisfy 

Pn,d = Xn ■ Pn,d-1 + Pn-l,d- (10) 

Along with the equations P^^i = clqXq + aixi + ■ ■ ■ + a^Xk and Po,<5 = PqXq, this gives a 
formula for the polynomial Pn^d- Clearly, some P^^s may be the zero polynomial. 

The rest of the proof is devoted to compute a bound on the size of the formula obtained 
by Equation (10). Let F{n,d) denote the bound on the size of the formula computing 
Pn,d: F{n,d) < F{n -l,d) + F{n,d- 1) + 2. For the base cases, F{k, 1) < k for aU k, 
f(o, 5) <S-1. Let G{N, d) = F{N - d - l,d) + 2 {for N > d and d > 1). Then G{N, d) 
satisfies Pascal's formula 

G{N,d)<G{N -l,d) + G{N -l,d-l) (11) 
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and G(<5 + 1, (5) < (5 + 1, G{k + 2,1) <k + 2. Thus G{N, d) is bounded from above by the 
binomial coefficient ( , ) , so we obtain 

F(n,d)<("+;;+^)-2. (12) 

This gives a first bound on F[n, d), somewhat bigger than the one announced. This comes 
from the fact that the base case bound F{0, 5) < 5 — 1 is too large: As the new variable 
xq is for homogenization, the actual formula is obtained by replacing it by 1 and therefore 
the formula for Pqs = Po^o is made of a single input labelled by 1 with the constant po on 
the wire going from it. So F{0, 6) = 0. 

This remark yields the same equation as Equation (11) for G but with a new base case 
G{6 + 1,5) = 2. A general form for such recurrences is 



G{N,d) = J2-.^/_^) 



3=0 

for some aj. Nevertheless, the values we get for the Oj if we apply this equation to the base 
cases are not really explicit. Therefore, we shall proceed in a different way: the new bound 
for G{N,d) is computed as the difference between the bigger bound (^) and the number 
of Po,(5 that were counted. In the recurrence (10), consider the recursion tree: Suppose 
that the vertex corresponding to Pn,d-i is the left child of the vertex corresponding to Pn,d^ 
and Pn-i^d its right child. The root of the recursion tree corresponds to the output of the 
formula, and its leaves to some Pk^i or some Po,(5- The quantity to count is the number 
of leaves corresponding to some -Po,<5- ^ path from the root Pn,d to Po,(5 has to decrease 
the first argument from n to and the second from d to 5. In the recursion tree, this 
corresponds to a path going n times to the right and [d — 8) times to the left. Moreover, 
such a path finishes by a move from Pi^s to its right child i-b,<5j as Po,<5+i has no child. Let 
us define the set of strings Wi,j as 

Wi^j = {w ^ {L, R}* : \w\fi = i and \w\l = j} ■ 

The cardinality of Wij is {^) as an element of this set is determined by the i places for the 
letters R in a length-(i + j) word. As the path from Pn^d to Po,5 finishes by a right move, 
the number of Pqs occurring in the recursion tree is equal to the cardinality of Wn-i^d-S^ 
that is ("^„I^~ ). And for each Pq^s, the original bound counted {6 — 1) operations instead 
of zero. Thus, to get a tighter bound we have to subtract 



(«-i) „:r =»-i-') 



5=1 -^ -^ ' j=0 -^ ^ 



Let Mon;^ (resp. Mon--') be the set of all monomials in n variables of degree j (resp. at 
most j). Then Mon:^ has cardinality { ]~ ), and {d — j — 1)(" j~ ) is the cardinality of 
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the set {x^ Mon;^ '■0<p<d — j — 2} where x is a fresh variable. Thus, the sum over j of 

those quantities is the cardinahty of Mon-,^ , that is (" , 7 ) • This gives the first part of 

the theorem: 

, „ /'n + d+l\ fn + d- 1\ 
F{n,d) <( _ ) - f _ 1-2. 



n+ 1 



n+ 1 



In the rest of the proof, we shall give a bound on the dimensions of the matrix obtained by 
our construction of Section 2. 

In [Quarez 2008], the symmetric matrix that is built contains linear functions as entries 
(and not only variables and constants). Therefore, we now give a bound in that case to 
permit a tighter comparison between both methods. In other words, we suppose that the 
inputs of the formula are not only constants and variables, but also linear functions. This 
amounts to defining the size of the arithmetic formula uqXq + aiXi + • • • + a^x^ as instead 
of k. As in the previous paragraph, a direct computation where the bounds on the base 
cases are changed can be done but yields non explicit formulas. Therefore, we use the same 
technique as before: The size of the formula when inputs can be linear functions is the 
difference between the size of the classical formula and the number of linear functions that 
appear. Those linear functions are the P^ i and appear as leaves in the recursion tree. A 
leaf labelled by Pk^i is reachable by a path going (n — k) times to the right and {d— 1) times 
to the left. As above, the path finishes by a move from P^ 2 to its left child Pfc,i- Therefore 
the number of leaves labelled by P^^i is the cardinality of Wn-k,d-2, that is (" n-k~ )■ ^^^ 
those leaves count for k additions, thus the total number of saved additions is 



k=l 



n + d — k 
n — k 



n-l 

j=0 



j + d- 2 



The computation is now the same as above and this sum equals ("^ ) 
Theorem 3, we get a symmetric matrix of dimensions 

^n + d+l\ /n + d—l\ fn + d—1 



S{n,d) <2 



n+ 1 



n+l 



n 



1 



To complete the proof, it is sufficient to use Pascal's formula twice: 



n + d+l 
n+l 



n + d\ fn + d 
n+l/ \ n 



n + d-l 
n+l 



+ 



n + d-l 



n 



+ 



n / \ n + 1 



n + d-l 
n 
n + d-l 
n — 1 



Using now 



n + d-l 
n-l 



U 



30 



Note that the bound F{n, d) we obtain with this construction is only better by a hnear 
factor in n than the obvious formula consisting of a sum of all the monomials. Indeed, for 
any J < d, there are at most (" v~ ) monomials of degree j which use (j — 1) multiplications, 

and there are at most (("J )"-'-) additions. Therefore the size of the formula we get in 
this way is 

E, ■ -, N /" + J ~ 1\ fn + d\ ^ fn + d\ nin + d) (n + d-l 

The first equality comes from similar techniques as in the previous proof and the second 
one is a straightforward computation. This yields a matrix of dimensions ^^ A S{n,d) 
approximately. 

Nevertheless, this is a bound in the worst case, that is for a polynomial M„^^ in which 
all the monomials of degree at most d appear. But in this special case one can change this 
construction if the aim is to have the polynomial M^^d itself. Indeed, the recurrence given 
by Equation (10) can be altered in the following manner: 

Mn,d = XnMn,d-l + Mn-l,d 

= XnMn4-l + Xn-lMn-14-l + Mn-2,d 
= XnMn4-l -\ h XoMo,d_l. 

This gives an inductive construction of a skew circuit to compute M^^d- At step 1, M„ 1 
is built, and it is clear that every M^-k,! is represented by a gate in the circuit. At step 
6 < d, suppose that we have a circuit such that every Mn-k,5-i is represented by a gate. 
Then one can build a circuit with (n + 1) new variable inputs, (n + 1) multiplication gates 
and n addition gates such that every Mn-k,5 is represented by a gate. At each step, the 
circuit size increases by (2n + 1) and (n + 1) inputs are added. As the size of the circuit for 
degree 1 is n with (n + 1) inputs, the circuit for Mn^d has size {2nd — n + d — 1) and has 
(n + l)d inputs. This yields a matrix of polynomial dimensions (in n and d), much smaller 
than with Quarez's construction. 

Let us now compare the bounds of Theorem 7 in the worst case with Quarez's. To 
this end let us consider a polynomial with n variables and of degree 2d. Then Quarez 
builds a symmetric matrix of dimensions 2 ("^ ) whereas our construction yields a matrix 
of dimensions 4("'^~ ) ~ ^- ^ bound on the quotient of those quantities can be given 
using the inequalities (see e.g. [Knuth 1997]) 

n + d\ fe{n + d)Y a f"'^'^^~^\ (n + 2d-l 



n J \ n J \'^/\ "^ 

So, the quotient is bounded by 

'e(n + (i)Y" / n \"^ _ n f n + d 
~n 7 ' \n + 2d-l) ~^ ' \n + 2d-l 



<e". 
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This means that Quarez's construction is exponentiahy better in the general case even 
though our construction yields much smaller matrices when the polynomial has a polynomial 
size formula or weakly-skew circuit. 

We now compare Quarez's results and ours for the special case of the permanent. This is 
an important example of a polynomial for which no polynomial size circuit is known (even 
non weakly-skew). Nevertheless, there exist formulas for computing it of much smaller 
size than the bounds for the general case [Ryser 1963; Glynn 2010]. For instance, Ryser's 
formula to compute the permanent of a matrix M is 

n 

per(A)= Y. (-l)""nE^^^- 
5C{l,...,n} i='^j<^S 

As the sums of variables are not counted, this gives a size-0(n2") formula, and hence yields 
a symmetric matrix of dimensions 0(n2") to represent the permanent. Let us consider the 
permanent of a (2n x 2n) matrix. This is a polynomial of degree 2n with An? variables. 
Therefore, Quarez's construction yields a matrix of dimensions 2( "„ ) • This quantity can 
be bounded as above and therefore we get the following bound (up to a constant factor) 
for the quotient: 



V")^ ((4n^+n)/n)" ^ 4-n" _^^, 



n22" - n4" ~ n4" 

A more careful computation via Stirling's formula shows that this quotient is equal to 
0(n"'~^"(4e)") when n tends to infinity. 

5 Characteristic 2 

In characteristic 2, the constructions of Sections 2 and 3 are not valid anymore because of 
the coefficients 1/2 they use. Nevertheless, for a polynomial computable by a weakly-skew 
circuit, it is possible to represent its square as the determinant of a symmetric matrix. On 
the other hand, representing the polynomial itself seems to be a challenging problem. For 
instance, it is not possible to represent the polynomial xy+z this way (Grenet et al. [2011b]), 
but we don't have for the moment a characterisation of representable polynomials. Related 
to these problems, the VNP-completeness of the partial permanent is also studied. Actually, 
we give an almost complete answer to an open question of Biirgisser [2000] (Problem 3.1) 
showing that if the partial permanent is complete in finite fields of characteristic 2, then the 
(boolean) polynomial hierarchy collapses. For any field of characteristic 2 (finite or infinite), 
we show that the VNP-completeness of this family would imply that every VNP family of 
polynomials has its square in VP^^ (i.e. has polynomial size weakly-skew circuits). This 
also seems unlikely to happen unless VP^^ = VNP. We refer to [Biirgisser 2000; Malod and 
Portier 2008] for the formal definitions of the complexity classes VNP and VPws- 

Let G be an edge- weighted graph with vertices {vi, . . . ,Vn}- Recall that the adjacency 
matrix ^ of G is the (n x n) symmetric matrix defined by Aij = Aji = Wij where Wij 
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is the weight of the edge ViVj. Suppose now that G is bipartite with two independent 
sets of vertices Vr and Vc of cardinahty m and n respectively. Let Vr = {ri, . . . ,rm} and 
Vc = {ci, . . . ,Cn}- The biadjacency matrix of G (also known as the bipartite adjacency 
matrix) is the (m x n) matrix B such that Bij is the weight of the edge between Tj and Cj. 
This means that the rows of B are indexed by Vr and its columns by Vc- For a bipartite 
graph G of adjacency and biadjacency matrices A and B respectively, 

B 
S* 

Throughout this section, we shall use the usual definition of the weight of a partial 
matching: it is the product of the weights of the edges it uses. 

5.1 Symmetric determinantal representation of the square of a polyno- 
mial 

Lemma 6. Let G be an edge-weighted graph and A its adjacency matrix. In characteristic 
2, the determinant of A is the sum of the weights of the cycle covers with cycles of length 
at most 2. 

Proof. Let us consider G as a symmetric digraph (that is an edge uv is seen as both arcs 
(u, v) and {v,u)). In Lemma 1, the signs of the cycle covers are considered. In characteristic 
2, this is irrelevant. Therefore, the determinant of A is the sum of the weights of the cycle 
covers of G. 

Let C be a cycle cover of G containing a (directed) cycle of length at least 3 denoted by 
{vi,V2, ■ ■ ■ ,Vk,vi). One can change the direction of this cycle (as G is symmetric) and obtain 
a new cycle cover G' containing the same cycles as C, but (vk,Vk-i, ■ ■ ■ ,vi,Vk) instead of 
(t"!, 112) ■ ■ ■ ,Vk,vi). Clearly, the weights of G and C" are the same as the graph is symmetric. 
Therefore, when the determinant of A is computed in characteristic 2, the contributions of 
those two cycle covers to the sum cancel out. This shows that the determinant of a matrix 
in characteristic two is obtained as the sum of the weights of cycle covers with cycles of 
length 1 (loops) or 2. D 

Proposition 1. Let p be a polynomial over a field of characteristic 2, represented by a 
weakly-skew circuit of fat size m. Then there exists a symmetric matrix A of dimensions 
(2m + 2) such that p^ = det(A). 

Proof. Let C be a weakly-skew circuit representing a polynomial p over a field of charac- 
teristic 2. Let M be the matrix obtained by Malod and Portier's construction [2008] such 
that perM = p. Let G be the digraph represented by M, and let G' be the bipartite graph 
obtained from G by the two following operations: Each vertex u of G is turned into two 
vertices v^ and v* in G', and each arc {u.,v) is turned into the edge {;u^,v^}. A loop on 
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a vertex u is simply represented as the edge {ii'*,n*}. Let A be the symmetric adjacency 
matrix of G' (when the vertices are ordered VQ,vf, . . . , v^, Vq, . . . , v^). 

It is weh-known that cycle covers of G and perfect matchings of G' are in one-to-one 
correspondence. If there is a cycle cover of G, then each vertex v belongs to a cycle, and 
thus has both a predecessor v and a successor w. This means that u* and u^ are matched to 
v^ and w* respectively (if u is covered by a loop, then u^ and u* are matched). Conversely, 
suppose that G' has a perfect matching. Let u'^ be any vertex. Then it is matched to some 
u*. In the same way, v^ is matched to some w^. As the set of vertices is finite, as some 
point we go back to u*. Thus it defines a cycle in G, and by doing the same process with 
other vertices not in this cycle this eventually defines a cycle cover in G. 

This one-to-one correspondence shows that the determinant of M equals the sum of the 
weights of the perfect matchings in G' . If a perfect matching in G' is considered as a cycle 
cover with length-2 cycles, the weight of the cycle cover is the square of the weight of the 
perfect matching. Indeed, in the cycle cover, all the arcs of the length-2 cycles have to be 
considered, that is each edge contributes twice to the product. Lemma 6 and the fact that 
there is no loop in G' show that 



det{A) = J2H^^? = {Y.''i^^^ 



where n ranges over all perfect matchings of G' and w{fi) is the weight of the perfect 
matching fi. The second equality holds as the field has characteristic 2. 

Finally, it is shown in [Malod and Portier 2008] that p = det(M), and we showed that 
det(M) = Em^(/^) and det(^) = (E/.^(m))^- Therefore, det(^) = det(M)2 =p'^. D 

This proposition raises the following question: Let / be a family of polynomials such 
that /^ E VP^s. Does / belong to VP^^? This question is discussed with more details in 
the next section. 

5.2 Is the partial permanent complete in characteristic 2? 

Definition 8. Let X = (Xij) be an (n x n) matrix. The partial permanent of X, as defined 
by Biirgisscr [2000], is 

per*(X) = j; H X,,(,), 

■T jgdcf(7r) 

where the sum ranges over the injective partial maps from [n] = {1, . . . , n} to [n] and def (vr) 
is the domain of the partial map vr (recall that a partial map is a map from a subset of [n] 
to [n]). 

The family (PER* ) is the family of polynomials such that PER* is the partial permanent 
of the (n X n) matrix whose coefficients are the indeterminates Xij . 
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Lemma 7. Let G be the complete bipartite graph with two independent sets of vertices Vr 
and Vc such that the edge between ri and Cj is labelled by Bij (the matrix B is the biadjacency 
matrix of G). Then the partial permanent of B is equal to the sum of the weights of the 
partial matchings of G. 

A partial matching in a graph G is a set of pairs of vertices connected by an edge such 
that no vertex appears in more than a pair. Equivalently, a partial matching can be seen 
as a set of disjoint edges. The weight of a partial matching is the product of the weights of 
its edges. 

The proof of the lemma is quite straightforward as a partial injective map vr from [n] 
to [n] exactly defines a partial matching in G such that for i € def(7r), rj is matched with 

Lemma 8. Let G be the complete bipartite graph with two independent sets of vertices Vr 
and Vc such that the edge between ri and Cj is labelled by Bij (the matrix B is the biadjacency 
matrix of G). Let A he its adjacency matrix. Then in characteristic 2, 

det(A + /2n) = (per*(5))^ 

where l2n is the identity matrix of dimensions 2n. 

Proof. By Lemma 6, to compute a determinant in characteristic 2, one can focus only on 
cycles of length at most 2. A cycle cover with such cycles actually is a partial matching 
when the graph is symmetric (length-2 cycles define the pairs of vertices, and length- 1 cycles 
are isolated vertices). Considering G as a symmetric digraph, the weight of a cycle cover is 
equal to the product of the weights of its loops and the square of the weights of the edges 
it uses (a length-2 cycle corresponds to an edge). 

Consider the graph G' obtained from G by adding weight-1 loops on all its vertices. In 
other words, G' is the graph whose adjacency matrix is ^ + /2„. By the previous remark, 
and by the fact that the loops have weight 1 , the determinant of ^4 + /2n is 

det(^ + L2n) = Y^ w{fif = (J2 w{l^) 

where /x ranges over the partial matchings of G' and u)(//) is the weight of the partial 
matching /i. The second equality is true as the characteristic of the field is 2. 

Recall now that G is bipartite. Of course, the partial matchings of G and G' are the 
same. So 

per* (5) =^u;(/i), 

where fi ranges over the partial matchings of G. This proves the lemma. D 

35 



An alternative proof of this lemma was suggested by an anonymous referee. In any 
field, we have the polynomial identity per(yl + tl2n) = X]fc=o''fc^^"~'^' '^^ere c^ is the sum 
of the permanents of all central minors of A of size k. In particular, we have det(A + /2n) = 
X]fc=o ''*: ^^ characteristic 2. A nonzero permanent of a central minor of A is of the form 



per 



■ M' 
M^ 



per(M)^ 



where M is a square submatrix of B. Hence det(^ + l2n) = X^a/cb P^^(-^)^' where M 
ranges over all square submatrices of B. Since we are in characteristic 2, this last sum is 

equal to (X^a/cb P^^(-^))^- -^^^ '^^ have per*(i?) = X]a/c_b P^'^C-^^) '^y definition of the 
partial permanent, and the conclusion of the lemma follows. 

Lemma 8 shows in particular that to compute the parity of the number of partial 
matchings in a bipartite graph, it is sufficient to compute a determinant (this is the case 
where G is not edge- weighted) . Therefore, this problem is solvable in polynomial time. 
This was already mentioned by Valiant [2005] but without any proof or reference. 

Theorem 8. In characteristic 2, the family ((PER*)^) is in VP^^. 

Proof. The previous lemma shows that the polynomial (PER*)^ is a p-projection of DET2n 
in characteristic 2. Thus, ((PER*)^) is in VP^^. D 

Suppose that (PER*) is VNP-complete. Then every VNP family (/„) is a p-projection 
of (per;), and thus (/^j is a p-projection of ((PER*)^). Let VNP^ = {(/2) : (/„) g VNP} 
be the class of squares of VNP families. This implies the following corollary of the theorem: 

Corollary 1. In any field of characteristic 2, if (PER^) is VHP -complete, then VNP C 

This situation is unlikely to happen. In particular, it would be interesting to investigate 
whether this inclusion implies that VP^^ = VNP in characteristic 2. Let us now give 
another consequence of (PERJj) being VNP-complete. This only holds for ffiiite fields of 
characteristic 2 but may give a stronger evidence that (PER.^) is unlikely to be VNP- 
complete. 

Theorem 9. // the partial permanent family is VHP -complete in a finite field of charac- 
teristic 2, then ©P/poly = NC /poly, and the polynomial hierarchy collapses to the second 
level. 

The proof of this theorem uses the boolean parts of Valiant's complexity classes defined 
in [Biirgisser 2000]. In the context of finite fields of characteristic 2, the boolean part 
of a family (/„) of polynomials with coefficients in the ground field F2 is the function 
bpf : {0, 1}* -^ {0, 1} such that for x G {0, 1}", bpf{x) = fn{x) (mod 2). The boolean part 
BP(C) of a Valiant's class C is the set of boolean parts of all f ^ C. 
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Proof. Let (/„) be a VNP family and {(fn) its boolean part. As ipnix) € {0,1} for all 
X G {0, 1}", (ipn) is the boolean part of (/2) too. This shows that BP(VNP) C BP(VNp2). 
By Corollary 1, VNP^ C VP^^ C VP. Thus, BP(VNP) C BP(VNp2) C BP(VP) and as 
VPC VNP 

BP(VP) = BP(VNP). 

Biirgisser [2000] shows that in a finite field of characteristic 2, ©P/poly = BP(VNP), 
and BP(VP) C NCVpoly. Hence, eP/poly C NCVpoly. Moreover, NCVpoly C P/poly C 
0P/poly hence we conclude that 

eP/poly = NCVpoly. 

The collapse of the polynomial hierarchy follows from a non uniform version of the 
Valiant- Vazirani Theorem [1986]: Theorem 4.10 in [Biirgisser 2000] states that N P/poly C 
®P/poly. Therefore, 

NCVpoly C NP/poly C eP/poly = NCVpoly. 

In particular, P/poly = NP/poly and Karp and Lipton [1982] showed that this implies the 
collapse of the polynomial hierarchy to the second level. D 

Since the submission of this paper, Biirgisser's open problem has been completely 
settled. Guillaume Malod [2011] has proved, using clow sequences a la Mahajan and 
Vinay [1997], that PER* S VP^^- Stefan Mengel subsequently noticed that the result 
can be derived from a result of Valiant on Pfaffian Sums ]Valiant 2002], see also [Guo et al. 
2011]. 

6 Conclusion 

Figure 18 shows the graphs obtained from the weakly-skew circuit and the formula of 
Fig. 1(b) and (c) for a field of characteristic different from 2, and Table 2 recalls all the 
constructions used in this paper. 

Table 1 compares the results obtained, in this paper and in previous ones. The bounds 
are given for a formula of green size e and for a weakly-skew circuit of green size e with i 
input gates labelled by a variable. 
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Figure 18: Graphs obtained from the weakly-skew circuit and the formula given in Fig. 1(b) 
and fc). 





Non-symmetric 
matrix 


Symmetric 
matrix 


Formula 


e + 1 


2e-hl" 


Weakly-skew circuit 


{e + i) + l 


2(e + i) + l 



"The bound is achieved if and only if the entries can be complex numbers. Else, the bound is 2e -I- 2. 

Table 1: Bounds for determinantal representations of formulas and weakly-skew circuits. 
The bounds for symmetric representations are new, and the bound for a non-symmetric 
representation of a weakly-skew circuit is a slight improvement of known bounds. 



38 



Valiant's construction 

with constants 

(Section 2.1) 



Formulas 



(Section 2.2) 



Formulas 

with constants 

(Section 2.2) 



Weakly skew 

circuits 
(Section 3.1) 



Weakly skew circuits 

with constants 

(Section 3.3) 






X ^ 



Input gate 



H 
constant 1 



X 

no constant 



X 

H 

constant 1 



f s 



• ^c 



no constant 



X 



G ] 

f V y 



-1 



-1 



; G' ] 
constants 1 



00 

CO 



Addition gate 






1 a 1 



I G' 



C2/C1 



tR\ 



V, 




1 






2' 
t-1 



no constant 




I — 1 



constants 1 



Multiplication 
gate 




2 1 



G2 ] 



constant C1C2 



no constant 



constant C1C2 



ua — tp 

no constant 




tp 



constant ciC2Cpc^ 



Table 2: Summary of the constructions 



The (e + 1) bound for the representation of a formula by a (non-symmetric) matrix 
determinant was given in [Liu and Regan 2006] by a method purely based on matrices. We 
show in Section 2.1 that this bound can also be obtained directly from Valiant's original 
proof when we remove the little flaw it contains. The (e + i + l) bound for the representation 
of a polynomial computed by a weakly-skew circuit can be obtained from the (m+l) bound 
(where m is the fat size of the circuit) obtained in [Malod and Portier 2008] if we use our 
minimization lemma (Lemma 5) as well as a similar trick as in the proof of Theorem 6. 
Both bounds for the symmetric cases are given in this paper. 

A formula is a special case of weakly-skew circuit. If our construction for weakly-skew 
circuits is applied to a formula, this yields a matrix that can be as large as twice the size of 
the matrix obtained with the specific constructions for the formulas. In the converse way, 
one could turn a weakly-skew circuit into a formula and then apply the construction for the 
formula. Yet, turning a weakly-skew circuit into a formula of polynomial size is not known 
to be possible. In fact, this would give a polynomial size formula for the determinant, and 
hence a parallel time upper bound of O(logn). So far, the best upper bound is Csansky's 
famous 0(log n) upper bound [1976]. 

All of these results are valid for any field of characteristic different from 2. We showed 
that there are some important differences for the complexity of polynomials over fields of 
characteristic 2. The question of characterizing which polynomials can be represented as 
determinants of symmetric matrices is quite intriguing and remains open. 
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